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Abstract. We describe the implementation of a subfield of the held of formal Puiseux 
series in polymake. This is employed for solving linear programs and computing convex 
hulls depending on a real parameter. Moreover, this approach is also useful for computations 
in tropical geometry. 


1 Introduction 

It is well known and not difficult to see that the standard concepts from linear programming 
(LP), e.g., the Farkas Lemma and LP duality, carry over to an arbitrary ordered field; e.g., see [7, 
Section II] or [16, §2.1]. Traces of this can already be found in Dantzig’s monograph ]8, Chapter 
22]. This entails that any algorithm whose correctness rests on these LP corner stones is valid 
over any ordered field. In particular, this holds for the simplex method and usual convex hull 
algorithms. A classical construction, due to Hilbert, turns a field of rational functions, e.g., with 
real coefficients, into an ordered field; see ]30, §147]. In ]16] Jeroslow discussed these fields in the 
context of linear programming in order to provide a rigorous foundation of the so-called “big M 
method”. The purpose of this note is to describe the implementation of the simplex method and of 
a convex hull algorithm over fields of this kind in the open source software system polymake ]14]. 

Hilbert’s ordered field of rational functions is a subfield of the field of formal Puiseux series 
with real coefficients. The latter field is real-closed by the Artin-Schreier Theorem [27, 
Theorem 12.10]; by Tarski’s Principle (cf. [28]) this implies that has the same first order 

properties as the reals. The study of polyhedra over motivated by tropical geometry 

]9], especially tropical linear programming ]2]. The connection of the latter with classical linear 
programming has recently lead to a counter-example ]1] to a “continuous analogue of the Hirsch 
conjecture” by Deza, Terlaky and Zinchenko ]10]. In terms of parameterized linear optimization 
(and similarly for the convex hull computations) our approach amounts to computing with 
sufficiently large (or, dually, sufficiently small) positive real numbers. Here we do not consider 
the more general algorithmic problem of stratifying the parameter space to describe all optimal 
solutions of a linear program for all choices of parameters; see, e.g., [17] for work into that 
direction. 

This paper is organized as follows. We start out with summarizing known facts on ordered 
fields. Then we describe a specific field, Q{t}, which is the field of rational functions with rational 
coefficients and rational exponents. This is a subfield of which we call the field of Puiseux 

fractions. It is our opinion that this is a subfield of the formal Puiseux series which is particularly 
well suited for exact computations with (some) Puiseux series; see [22] for an entirely different 
approach. In the context of tropical geometry Markwig ]23] constructed a much larger field, which 
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contains the classical Puiseux series as a proper subfield. For our applications it is relevant to 
study the evaluation of Puiseux fractions at sufficiently large rational numbers. In Section 3 
we develop what this yields for comparing convex polyhedra over with ordinary convex 

polyhedra over the reals. The tropical geometry point of view enters the picture in Section 4. 
We give an algorithm for solving the dual tropical convex hull problem, i.e., the computation of 
generators of a tropical cone from an exterior description. Allamigeon, Gaubert and Goubault 
gave a combinatorial algorithm for this in [4], while we use a classical (dual) convex hull algorithm 
and apply the valuation map. The benefit of our approach is more geometric than in terms of 
computational complexity: in this way we will be able to study the fibers of the tropicalization 
map for classical versus tropical cones for specific examples. Section 5 sketches the polymake 
implementation of the Puiseux fraction arithmetic and the LP and convex hull algorithms. The 
LP solver is a dual simplex algorithm with steepest edge pivoting, and the convex hull algorithm 
is the classical beneath-and-beyond method [11] [18]. An overview with computational results is 
given in Section 6. 

Acknowledgment We thank Thomas Opfer for contributing to and maintaining within the 
polymake project his implementation of the dual simplex method, originally written for his 
Master’s Thesis [25]. 

2 Ordered fields and rational functions 

A field F is ordered if there is a total ordering < on the set F such that for all a,b,c G F the 
following conditions hold: 

(i) if a < & then a + c < b + c, 

(ii) if 0 < a and 0 < b then 0 < a ■ b. 

Any ordered field necessarily has characteristic zero. Examples include the rational numbers Q, 
the reals K and any subfield in between. 

Given an ordered field F we can look at the ring of univariate polynomials F[t] and its quotient 
field F(t), the field of rational functions in the indeterminate t with coefficients in F. On the ring 
F[t] we obtain a total ordering by declaring p < q whenever the leading coefficient of g — p is a 
positive element in F. Extending this ordering to the quotient field by letting 

u p 

— < - : uq <vp , 

V q 

where the denominators v and q are assumed positive, turns F(t) into an ordered field; see, e.g., 
[30, §147]. This ordered field is called the “Hilbert field” by Jeroslow [16[. 

By definition, the exponents of the polynomials in F[t] are natural numbers. However, 
conceptually, there is no harm in also taking negative integers or even arbitrary rational numbers 
as exponents into account, as this can be reduced to the former by clearing denominators and 
subsequent substitution. For example, 

2t3/2_t-i _ 2^5/2 _i _ 2 s15-1 

1 + 3t-i/3 “ t + 3f2/3 “ s6 + 3 s4 ’ ^ ^ 

where s = In this way that fraction is written as an element in the field of rational 

functions in the indeterminate s = with rational coefficients. Further, if p S F(t^/“) and 
q G for natural numbers a and (3, then the sum p + q and the product p • q are contained 

in F(f^/scd(“./3)This shows that the union 

F{t} = |jF(ti/^) 

i/>i 


( 2 ) 
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is again an ordered field. We call its elements Puiseux fractions. The field F{t} is a subfield of 
the field of formal Puiseux series, i.e., the formal power series with rational exponents of 

common denominator. For an algorithmic approach to general Puiseux series see [22]. 

The map val which sends the rational function p/q, where p,q G F[t^/'^], to the number 
degtP — degt q defines a non-Archimedean valuation on F(t). Here we let val(O) = oo. As usual 
the degree is the largest occurring exponent. The valuation map extends to Puiseux series. More 
precisely, for f,g G F{t} we have the following: 

(i) val(/ • g) = val(/) -fi val]^), 

(ii) val(/ + g) < max(val(/), val]^)). 

If F = K is the field of real numbers we can evaluate a Puiseux fraction / G K{t} at a real 
number r to obtain the real number /(r). This map is defined for all r > 0 except for the finitely 
many poles, i.e., zeros of the denominator. Restricting the evaluation to positive numbers is 
necessary since we are allowing rational exponents. The valuation map satisfies the equation 

hm logJ/(r)| = val(/) . (3) 

T—^OO 

That is, seen on a logarithmic scale, taking the valuation of / corresponds to interpreting t like 
an infinitesimally large number. Reading the valuation map in terms of the limit (3) is known as 
Maslov dequantization, see [24]. 

Occasionally, it is also useful to be able to interpret t as a small infinitesimal. To this end, 
one can define the dual degree deg*, which is the smallest occurring exponent. This gives rise to 
the dual valuation map Yal*{p/q) = degf p — deg^ q which yields 

val*(/-fi 5 ) > min(vaT(/),var(g)) and lim log^ |/(t)| = vaT(/) . 

T—>0 

Changing from the primal to the dual valuation is tantamount to substituting t by t~^. 

Remark 1. The valuation theory literature often employs the dual definition of a valuation. The 
equation (3) is the reason why we usually prefer to work with the primal. 

Up to isomorphism of valuated fields the valuation on the field ¥{t) of rational functions is 
unique, e.g., see [30, §147]. As a consequence the valuation on the slightly larger field of Puiseux 
fractions is unique, too. 

To close this section let us look at the algorithmically most relevant case F = Q. Then, in 
general, the evaluation map sends positive rationals to not necessarily rational numbers, again due 
to fractional exponents. By clearing denominators in the exponents one can see that evaluating 
at cr > 0 ends up in the totally real number field Q( -tya) for some positive integer v. For instance, 
evaluating the Puiseux fraction from Example (1) would give an element of Q(v^). 


3 Parameterized Polyhedra 

Consider a matrix A G Then the set 

C := {xGF{t}'^+i I A-a;>0} 

is a polyhedral cone in the vector space F{f}‘^+^. Equivalently, C is the set of feasible solutions of 
a linear program with d-|- 1 variables over the ordered field F{f} with m homogeneous constraints, 
the rows of A. The Farkas-Minkowski-Weyl Theorem establishes that each polyhedral cone is 
finitely generated. A proof for this result on polyhedral cones over the reals can be found in 
[31, §1.3 and §1.4] under the name ’’Main theorem for cones”. It is immediate to verify that the 
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arguments given hold over any ordered field. Therefore, there is a matrix B G for 

some n € N, such that 

C = {B-aj aeF{t}”, a > 0} . (4) 

The columns of B are points and the cone C is the non-negative linear span of those. 

Let L be the lineality space of C, i.e., L is the unique maximal linear subspace of 
which is contained in C. If dim L = 0 the cone C is pointed. Otherwise, the set C/L is a pointed 
polyhedral cone in the quotient space ¥{t}‘^~^^/L. A face of C is the intersection of C with a 
supporting hyperplane. The faces are partially ordered by inclusion. Each face contains the 
lineality space. Adding the entire cone C as an additional top element we obtain a lattice, the 
face lattice of C. The maximal proper faces are the facets which form the co-atoms in the face 
lattice. The combinatorial type of C is the isomorphism class of the face lattice (e.g., as a partially 
ordered set). Notice that our definition says that each cone is combinatorially equivalent to its 
quotient modulo its lineality space. 

Picking a positive element r yields matrices A{t) G and B{t) G as well 

as a polyhedral cone C{t) = {x G F'^+^ | A(t) • a: > 0} by evaluating the Puiseux fractions at the 
parameter t. Here and below we will assume that r avoids the at most finitely many poles of the 
(m -I- n) • (d -I- 1) coefficients of A and B. 

Theorem 2. There is a positive element tq S F so that for every t > tq we have 

C{t) = { B{t) • a I a G F”, a > 0} , 

and evaluating at t maps the lineality space of C to the lineality space of C (r). Moreover, the 
polyhedral cones C and C{t) over¥{t} and¥, respectively, share the same combinatorial type. 

Proof. First we show that an orthogonal basis of the lineality space L evaluates to an orthogonal 
basis of the lineality space of C(r). For this, consider two vectors x,y G F{t}‘^+^ and pick r large 
enough to avoid their poles and zeros. Then, the scalar product of x and y vanishes if and only if 
the scalar product of x{t) and yij) does. Hence, the claim follows. 

Now we can assume that the polyhedral cone C is pointed, i.e., it does not contain any linear 
subspace of positive dimension. If this is not the case the subsequent argument applies to the 
quotient CfL. 

Employing orthogonal bases, as for the lineality spaces above, shows that the evaluation maps 
the linear hull of C to the linear hull of C'(r), preserving the dimension. So we may assume that 
C is full-dimensional, as otherwise the arguments below hold in the linear hull of C. 

Let i < (™) be the number of d-element sets of linearly independent rows of the matrix A. 
For each such set of rows the set of solutions to the corresponding homogeneous system of linear 
equations is a one-dimensional subspace of For each such system of homogeneous linear 

equations pick two non-zero solutions, which are negatives of each other. We arrive at vectors 
in which we use to form the columns of the matrix Z G 

By the Farkas-Minkowski-Weyl theorem, we may assume that the columns of B from (4) only 
consist of the rays of C and that the rays of C form a subset of the columns of Z. In particular, 
the columns of B occur in Z. Since the cone C is pointed, the matrix B contains at most one 
vector from each opposite pair of the columns of Z. This entails that B has at most I columns. 

Further, the real matrix Z{t) contains all rays of C{t) for each r that avoids the poles of A 
and Z. In the following, we want to show that those columns of Z{t) which form the rays of C'(t) 
are exactly the columns of B{t). 

We define s{j, k) G F{t} to be the scalar product of the jth row of A and the fcth column of 
Z. The m • 2£ signs of the scalar products s(j, k), for j G [m\ and k G [2£], form the chirotope of 
the linear hyperplane arrangement defined by the rows of A (in fact, due to taking two solutions 
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for each homogenous system of linear equations, we duplicate the information of the chirotope). 
For almost all r S F evaluating the Puiseux fractions s{j, k) at r yields an element of F. For 
sufficiently large r the sign of s(j, k) agrees with its evaluation. This follows from the definition 
of the ordering on F{t}, cf. [16, Proposition, §1.3]. 

Let To G F be larger than all the at most finitely many poles of A and Z. Further, let tq be 
large enough such that the chirotope of A{t) agrees with the chirotope of A for all t > tq. 

By construction the rays of C correspond to the non-negative columns of the chirotope whose 
support, given by the 0 entries, is maximal with respect to inclusion; these are exactly the columns 
of B. The corresponding columns of the chirotope of A(t), for r > tq, yield the rays of C'(r), 
which, hence, are the columns of B{t). 

The same holds for the facets of C and C(r). The facets of C correspond to the non-negative 
rows of the chirotope whose support, given by the 0 entries, is maximal with respect to inclusion. 

Now the claim follows since the face lattice of a polyhedral cone is determined by the incidences 
between the facets and the rays. □ 

A statement related to Theorem 2 occurs in Benchimol’s PhD thesis [5]. The Proposition 5.12 
in [5] discusses the combinatorial structure of tropical polyhedra (arising as the feasible regions 
of tropical linear programs). Yet here we consider the relationship between the combinatorial 
structure of Puiseux polyhedra and their evaluations over the reals. As in the proof of [5, 
Proposition 5.12] we could derive an explicit upper bound on the optimal tq. To this end one 
can estimate the coefficients of the Puiseux fractions in Z, which are given by determinantal 
expressions arising from submatrices of A. Their poles and zeros are bounded by Cauchy bounds 
(e.g., see ]26, Thm. 8.1.3]) depending on those coefficients. We leave the details to the reader. 

A convex polyhedron is the intersection of finitely many linear inequalities. It is a called a 
polytope if it is bounded. Restricting to cones allows a simple description in terms of homoge¬ 
neous linear inequalities. Yet this encompasses arbitrary polytopes and polyhedra, as they can 
equivalently be studied through their homogenizations. In fact, all implementations in polymake 
are based on this principle. For further reading we refer to ]31, §1.5]. We visualize Theorem 2 
with a very simple example. 

Example 3. Consider the polytope P in for large t defined by the four inequalities 

Xi, X2 > 0, Xi + X2 <^1 Xi — X2 < t . 

The evaluations at r G {0,1,3} are depicted in Figure 1. For r = 0 we obtain a triangle, for r = 1 
a quadrangle and for r > 3 a triangle again. The latter is the combinatorial type of the polytope 
P over the field of Puiseux fractions with real coefficients. 

Corollary 4. The set of combinatorial types of polyhedral cones which can be realized over¥{t} 
is the same as over F. 

Proof. One inclusion is trivial since F is a subfield of F{t}. The other inclusion follows from the 
preceding result. □ 

For A G F{t}'"^‘^, b G Fit}"* and c G we consider the linear program LP(A, 6, c) over 

F{t} which reads as 

maximize ■ x 

subject to A-x = b,x>0. ^ 

For each positive r G F (which avoids the poles of the Puiseux fractions which arise as coefficients) 
we obtain a linear program LP(A(t), 6(t), c(r)) over F. Theorem 2 now has the following 
consequence for parametric linear programming. 
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Fig. 1. Polygon depending on a real parameter as defined in Example 3 


Corollary 5. Let x* G be an optimal solution to the LP (5) with optimal value v G F{t}. 

Then there is a positive element tq G F so that for every t > tq the vector x*(t) is an optimal 
solution for LP(A(r), 6(r), c(t)) with optimal value v(t). 

The above corollary was proved by Jeroslow [16, §2.3]. His argument, based on controlling 
signs of determinants, is essentially a local version of our Theorem 2. Moreover, determining all 
the rays of a polyhedral cone can be reduced to solving sufficiently many LPs. This could also be 
exploited to derive another proof of Theorem 2 from Corollary 5. 

Remark 6. It is worth to mention the special case of a linear program over the field F{t}, where 
the coordinates of the linear constraints, in fact, are elements of the field F of coefficients, but the 
coordinates of the linear objective function are arbitrary elements in F{f}. That is, the feasible 
domain is a polyhedron, P, over F. Evaluating the objective function at some r S F makes one 
of the vertices of P optimal. Solving for all values of r, in general, amounts to computing the 
entire normal fan of the polyhedron P. This is equivalent to solving the dual convex hull problem 
over F for the given inequality description of P; see also [17]. Here we restrict our attention to 
solving parametric linear programs via Corollary 5. 

The next example is a slight variation of a construction of Goldfarb and Sit [15]. This is a 
class of linear optimization problems on which certain versions of the simplex method perform 
poorly. 

Example 1. We fix d > 1 and pick a positive <5 < ^ as well as a positive e < |. Consider the 
linear program 

maximize 

subject to 0 < cci < 

Xj-i < Sxj < e‘^~^6 — Xj-i for 2 < j < d . 

The feasible region is combinatorially equivalent to the d-dimensional cube. Applying the simplex 
method with the “steepest edge” pivoting strategy to this linear program with the origin as the 
start vertex visits all the 2'^ vertices. Moreover, the vertex-edge graph with the orientation induced 
by the objective function is isomorphic to (the oriented vertex-edge graph of) the Klee-Minty 
cube [20]. See Figure 2 for a visualization of the 3-dimensional case. 

We may interpret this linear program over the reals or over (IR{i5}){£}, the field of Puiseux 
fractions in the indeterminate e with coefficients in the field M{(5}. This depends on whether we 
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Fig. 2. The 3-dimensional Goldfarb-Sit cube. 


want to view S and s as indeterminates or as real numbers. Here we consider the ordering induced 
by the dual valuation val*, i.e., S and e are small infinitesimals, where e S. Two more choices 
arise from considering e a constant in IR{(5} or, conversely, S a constant in K{e}. Note that our 
constraints on 6 and e are feasible in all four cases. 

Our third and last example is a class of linear programs occurring in [1]. For these the central 
path of the interior point method with a logarithmic barrier function has a total curvature which 
is exponential as a function of the dimension. 

Example 8. Given a positive integer r, we define a linear program over the field Q{t} (with the 
primal valuation) in the 2r -|- 2 variables uq, uq, ui, ui, . .., Ur, fr as follows: 

minimize Vg 
subject to Ug < t ^ Vg < 

Ui < tUi-i , Ui < tVi-i 

Vi < -I- 

Ur > 0 , Ur > 0 . 

Here it would be interesting to know the exact value for the optimal Tg in Theorem 2, as a 
function of r. Experimentally, based on the method described below, we found tq = 1 for r = 1 
and Tg = 2^ for r at most 5. We conjecture the latter to be the true bound in general. 

To find the optimal bound for a given constraint matrix A we can use the following method. 
One can solve the dual convex hull problem for the cone C, which is the feasible region in 
homogenized form, to obtain a matrix B whose columns are the rays of C. This also yields a 
submatrix of A corresponding to the rows which define facets of C. Without loss of generality we 
may assume that that submatrix is A itself. Let Tg be the largest zero or pole of any (Puiseux 
fraction) entry of the matrix A ■ B. Then for every value t > Tg the sign patterns of {A ■ B)(t) 
and A ■ B coincide, and so do the combinatorial types of C and C{t). Determining the zeros and 
poles of a Puiseux fraction amounts to factorizing univariate polynomials. 


for 1 < i < r 


4 Tropical Dual Convex Hulls 

Tropical geometry is the study of the piecewise linear images of algebraic varieties, defined over 
a field with a non-Archimedean valuation, under the valuation map; see [21] for an overview. 
The motivation for research in this area comes from at least two different directions. First, 
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tropical varieties still retain a lot of interesting information about their classical counterparts. 
Therefore, passing to the tropical limit opens up a path for combinatorial algorithms to be applied 
to topics in algebraic geometry. Second, the algebraic geometry perspective offers opportunities 
for optimization and computational geometry. Here we will discuss how classical convex hull 
algorithms over fields of Puiseux fractions can be applied to compute tropical convex hulls; see 
[19] for a survey on the subject; a standard algorithm is the tropical double description method 
of [3]. 

The tropical semiring T consists of the set K U {— 00 } together with u® v = max(u, v) as the 
addition and u Q v = u + v as the multiplication. Extending these operations to vectors turns 
T^+i g. semimodule. A tropical cone is the sub-semimodule 

tcone(G’) = { Ai © gi 0 ■ ■ • © I Ai,..., A^ G T} 

generated from the columns gi,... ,gn of the matrix G S Similar to classical cones, 

tropical cones admit an exterior description [13]. It is known that every tropical cone is the 
image of a classical cone under the valuation map val: —> T; see [9[. Based on this idea, we 

present an algorithm for computing generators of a tropical cone from a description in terms of 
tropical linear inequalities; see Algorithm 1 below. 

Before we can start to describe that algorithm we first need to discuss matters of general 
position in the tropical setting. The tropical determinant of a square matrix U € is given by 

tdet(f7) = 0 © ■ ■ ■ © . (6) 

o-eSe 


Here Si is the symmetric group of degree £', computing the tropical determinant is the same as 
solving a linear assignment optimization problem. Consider a pair of matrices H~ € 
which serve as an exterior description of the tropical cone 


Q 


z e 


r|p(d-|-l) 


H+ Qz> Qzj 


( 7 ) 


In contrast to the classical situation we have to take two matrices into account. This is due to 
the lack of an additive inverse operation. We will assume that p{i,j) ■= i 

any pair (f, j) € [m] x [d © 1], i.e., for each coordinate position at most one of the corresponding 
entries in the two matrices is finite. Then we can define 


Xihj) 


1 if p{i, j) = H± ^-00 
-1 if p{i, j) = Hr. ^-00 
0 otherwise . 


For each term Ui 7 r(i) © • ■ ■ © in (6) we define its sign as 

sgn(7r) •x(l,7r(l))---x(^,7r(^)) , 

where sgn(7r) is the sign of the permutation tt. Now the exterior description (7) of the tropical cone 
Q is tropically sign-generic if for each square submatrix U of iJ'*" © H~ we have tdet(f7) —00 
and, moreover, the signs of all terms Mi 7 r(i) © ■ • • © which attain the maximum in (6) agree. 

By looking at 1 x 1-submatrices U we see that in this case all coefficients of the matrix id + © H~ 
are finite and thus x(*ij) is never 0. 
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Algorithm 1: A dual tropical convex hull algorithm 

Input: pair of matrices H^,H~ G which provide a tropically sign-generic 

exterior description of the tropical cone Q from (7) 

Output: generators for Q 

pick two matrices A+, A~ G with strictly positive entries such that 

val(A+) = and val(A“) = H~ 

apply a classical dual convex hull algorithm to determine a matrix B G such 

that { i? • a I a G a > 0} = { x € | ~ ^~) • a; > 0, a; > O} 

return val(i?) 


Proof (Correctness of Algorithm 1). The main lemma of tropical linear programming [2, The¬ 
orem 16] says the following. In the tropically sign-generic case, an exterior description of a 
tropical cone can be obtained from an exterior description of a classical cone over Puiseux series 
by applying the valuation map to the constraint matrix coefficient-wise. This statement assumes 
that the classical cone is contained in the non-negative orthant. We infer that 


Q = i7+©z>iJ-©z} 

= vald^eKItl^^^-^+i) I A+-a: > A-> O}) 
= val({i?-a| a e ^ ^ O}) • 


Now [9, Proposition 2.1] yields Q = val({il • a | a S ^ 0}) = tcone(val(i3)). This ends 

the proof. □ 


The correctness of our algorithm is not guaranteed if the genericity condition is not satisfied. 
The crucial properties of the lifted matrices A+, A“ are not necessarily fulfilled. It is an open 
question of how an exterior description over T is related to an exterior description over in 

the general setting. We are even lacking a convincing concept for the “facets” of a general tropical 
cone. 


5 Implementation 

As a key feature the polymake system is designed as a Perl/C++ hybrid, that is, both pro¬ 
gramming languages are used in the implementation and also both programming languages can 
be employed by the user to write further code. One main advantage of Perl is the fact that 
it is interpreted; this makes it suitable as the main front end for the user. Further, Perl has 
its strengths in the manipulation of strings and file processing. C++ on the other hand is a 
compiled language with a powerful template mechanism which allows to write very abstract code 
which, nonetheless, is executed very fast. Our implementation, in C++, makes extensive use of 
these features. The implementation of the dual steepest edge simplex method, contributed by 
Thomas Opfer, and the beneath-beyond method for computing convex hulls (see [11] and [18]) 
are templated. Therefore polymake can handle both computations for arbitrary number field 
types which encode elements in an ordered field. 

Based on this mechanism we implemented the type RationalFunction which depends 
on two generic template types for coefficients and exponents. Note that the field of coefficients 
here does not have to be ordered. Our proof-of-concept implementation employs the classical 
Euclidean GCD algorithm for normalization. Currently the numerator and the denominator are 
chosen coprime such that the denominator is normalized with leading coefficient one. For the 
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most interesting case F = Q it is known that the coefficients of the intermediate polynomials 
can grow quite badly, e.g., see [12, Example 1]. Therefore, as expected, this is the bottleneck of 
our implementation. In a number field or in a field with a non-Archimedean valuation the most 
natural choice for a normalization is to pick the elements of the ring of integers as coefficients. The 
reason for our choice is that this more generic design does not make any assumption on the field of 
coefficients. This makes it very versatile, and it fits the overall programming style in polymake. 
A fast specialization to the rational coefficient case could be based on [12, Algorithm 11.4[. This 
is left for a future version. 

The polymake implementation of Puiseux fractions F{t} closely follows the construction de¬ 
scribed in Section 2. The new number type is derived from RationalFunction with overloaded 
comparison operators and new features such as evaluating and converting into TropicalNumber. 
An extra template parameter MinMax allows to choose whether the indeterminate t is a small or 
a large infinitesimal. 

There are other implementations of Puiseux series arithmetic, e.g., in Magma [6[ or MATLAB [29[. 
However, they seem to work with finite truncations of Puiseux series and floating-point coefficients. 
This does not allow for exact computations of the kind we are interested in. 

6 Computations 

We briefly show how our polymake implementation can be used. Further, we report on timings 
for our LP solver, tested on the Goldfarb-Sit cubes from Example 7, and for our (dual) convex 
hull code, tested on the polytopes with a “long and winded” central path from Example 8. 


6.1 Using polymake 

The following code defines a 3-dimensional Goldfarb-Sit cube over the field Q{t}, see Example 7. 
We use the parameters e = t and 6 = ^. The template parameter Min indicates that the ordering 
is induced by the dual valuation val*, and hence the indeterminate t plays the role of a small 
infinitesimal. 

polytope > $monomial=new UniMonomial<Rational, Rational> ( 1 ) ; 

polytope > $t=new PuiseuxFraction<Min> (Smonomial); 
polytope > $p=goldfarjb_sit (3,2*$t, 1/2) ; 

The polytope object, stored in the variable $p, is generated with a facet description from which 
further properties will be derived below. It is already equipped with a LinearProgram subobject 
encoding the objective function from Example 7. The following lines show the maximal value 
and corresponding vertex of this linear program as well as the vertices derived from the outer 
description. Below, we present timings for such calculations. 

polytope > print $p->LP->MAXIMAL_VALUE; 

( 1 ) 

polytope > print $p->LP->MAXIMAL_VERTEX; 


(1) 

(0) 

(0) (1) 




polytope > print 

$P- 

>VERTICES; 


(1) 

(0) 

(0) (0) 




(1) 

(t^' 

2) (2*t^2) 

(4 

*t^2) 


(1) 

(0) 

(t) (2*t) 




(1) 

(t^' 

2) (t -2*t 

-2) 

(2*t -4*t' 

'2) 

(1) 

(0) 

(0) (1) 




(1) 

(t- 

2) (2*t"2) 

(1 

-4*t^2) 


(1) 

(0) 

(t) (1 -2 

*t) 



(1) 

(t^' 

2) (t -2*t 

-2) 

(1 -2*t + 

4*t^2) 
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As an additional benefit of our implementation we get numerous other properties for free. For 
instance, we can compute the parameterized volume, which is a polynomial in t. 

polytope > print $p->VOLUME; 

(t"3 -4*t"4 + 4*t-5) 

That polynomial, as an element of the field of Puiseux fractions, has a valuation, and we can 
evaluate it at the rational number 

polytope > print $p->VOLUME->val; 

3 

polytope > print $p->VOLUME->evaluate(1/12); 

25/62208 

6.2 Linear programs 

We have tested our implementation by computing the linear program of Example 7 with polyhedra 
defined over Puiseux fractions. 

The simplex method in polymake is an implementation of a (dual) simplex with a (dual) 
steepest edge pricing. We set up the experiment to make sure our Goldfarb-Sit cube LPs behave 
as badly as possible. That is, we force our implementation to visit all n = 2'^ vertices, when d is 
the dimension of the input. Table 1 illustrates the expected exponential growth of the execution 
time of the linear program. In three of our four experiments we choose <5 as The computation 
over Q{a} costs a factor of about 80 in time, compared with the rational cubes for a modest e = g. 
However, taking a small £ whose binary encoding takes more than 18,000 bits is substantially 
more expensive than the computations over the field Q{e} of Puiseux fractions. Taking 5 as 
a second small infinitesimal is possible but prohibitively expensive for dimensions larger than 
twelve. 


Table 1. Timings (in seconds) for the Goldfarb-Sit cubes of dimension d with 5 = |. For e we tried a 
small infinitesimal as well as two rational numbers, one with a short binary encoding and another one 
whose encoding is fairly large. For comparison we also tried both parameters as indeterminates. 


d 

m 

n 

£ 

Q 

^=1 

Q 

2 

iQmie} 

^ ]^745UU 

3 

6 

8 

0.010 

0.003 

0.005 

0.101 

4 

8 

16 

0.026 

0.001 

0.017 

0.353 

5 

10 

32 

0.064 

0.002 

0.065 

1.034 

6 

12 

64 

0.157 

0.007 

0.253 

2.877 

7 

14 

128 

0.368 

0.006 

0.829 

7.588 

8 

16 

256 

0.843 

0.016 

2.643 

19.226 

9 

18 

512 

1.906 

0.039 

7.703 

47.806 

10 

20 

1024 

4.258 

0.090 

21.908 

118.106 

11 

22 

2048 

9.383 

0.191 

59.981 

287.249 

12 

24 

4096 

20.583 

0.418 

160.894 

687.052 


6.3 Convex hulls 

We have also tested our implementation by computing the vertices of the polytope from Example 8. 
For this we used the client long_and_winding which creates the d = (2r + 2)-dimensional 
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polytope given by to = 3r + 4 facet-defining inequalities. Over the rationals we evaluated the 
inequalities at 2^ which probably gives the correct combinatorics; see the discussion at the end of 
Example 8. This very choice forces the coordinates of the defining inequalities to be integral, such 
that the polytope is rational. The number of vertices n is derived from that rational polytope. 
The running times grow quite dramatically for the parametric input. This overhead could be 
reduced via a better implementation of the Puiseux fraction arithmetic. 


Table 2. Timings (in seconds) for convex hull computation of the feasibility set from Example 8. All 
timings represent an average over ten iterations. If any test exceeded a one hour time limit this and all 
larger instances of the experiment were skipped and marked —. 


r 

d 

m 

n 

Q{t} 

Q 

1 

4 

7 

11 

0.018 

0.000 

2 

6 

10 

28 

0.111 

0.000 

3 

8 

13 

71 

0.754 

0.010 

4 

10 

16 

182 

15.445 

0.036 

5 

12 

19 

471 

1603.051 

0.150 

6 

14 

22 

1226 

— 

0.737 

7 

16 

25 

3201 

— 

4.001 

8 

18 

28 

8370 

— 

25.093 

9 

20 

31 

21901 

— 

223.240 

10 

22 

34 

57324 

1891.133 


6.4 Experimental setup 

Everything was calculated on the same Linux machine with polymake perpetual beta version 2.15- 
betaS which includes the new number type, the templated simplex algorithm and the templated 
beneath-and-beyond convex hull algorithm. All timings were measured in CPU seconds and 
averaged over ten iterations. The simplex algorithm was set to use only one thread. 

All tests were done on openSUSE 13.1 (x86_64), with Linux kernel 3.11.10-25, clang 3.3 and 
perl 5.18.1. The rational numbers use a C++-wrapper around the GMP library version 5.1.2. As 
memory allocator polymake uses the pool_allocator from libstdc++, which was version 
4.8.1 for the experiments. 

The hardware for all tests was: 

Intel(R) Core(TM) i7-3930K CPU @ 3.20GHz 
bogomips: 6400.21 
MemTotal: 32928276 kB 
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